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ABSTRACT 



(N , 

Context. The source spectrum of cosmic rays is not well determined by diffusive shock acceleration models. The propagated fluxes 
' ' of proton, helium, and heavier primary cosmic-ray species (up to Fe) are a means to indirectly access it. But how robust are the 

' constraints, and how degenerate are the source and transport parameters? 

, Aims. We check the compatibility of the primary fluxes with the transport parameters derived from the B/C analysis, but also ask 

whether they add further constraints. We study whether the spectral shapes of these fluxes and their ratios are mostly driven by source 
or propagation effects. We then derive the source parameters (slope, abundance, and low-energy shape). 
r^' Methods. Simple analytical formulae are used to address the issue of degeneracies between source/transport parameters, and to 

Q , understand the shape of the p/He and C/0 to Fe/0 data. The full analysis relies on the USINE propagation package, the MINUIT 

5_j ■ minimisation routines (x^ analysis) and a Markov Chain Monte Carlo (MCMC) technique. 

Results. Proton data are well described in the simplest model defined by a power-law source spectrum and plain diffusion. They 
' can also be accommodated by models with, e.g., convection and/or reacceleration. There is no need for breaks in the source spectral 

indices below ~ 1 TeV/n. Fits to the primary fluxes alone do not provide physical constraints on the transport parameters. If we leave 
the source spectrum free, parametrised by the form dQ/dE = qP'^H^" , and fix the diffusion coefficient KiR) = Koj}''''^'R? so as to 
' reproduce the B/C ratio, the MCMC analysis constrains the source spectral index a to be in the range 2.2 - 2.5 for all primary species 

^ ' up to Fe, regardless of the value of the diffusion slope S. The values of the parameter 775 describing the low-energy shape of the source 

0^ , spectrum are degenerate with the parameter 777- describing the low-energy shape of the diffusion coefficient: we find rjs -tjt « for p 

00 . and He data, but 775 - ^/r ~ 1 for C to Fe primary species. This is consistent with the toy-model calculation in which the shape of the 

0^ ' P/He and C/0 to Fe/O data is reproduced if 77s - 777- » - 1 (no need for different slopes a). When plotted as a function of the kinetic 

' energy per nucleon, the low-energy p/He ratio is determined mostly by the modulation effect, whereas primary/0 ratios are mostly 

determined by their destruction rate. 

Conclusions. Models based on fitting B/C are compatible with primary fluxes. The different spectral indices for the propagated 
primary fluxes up to a few TeV/n can be naturally ascribed to transport effects only, implying universality of elemental source spectra. 



Key words. Methods: statistical - ISM: cosmic rays 



1. Introduction riety of possible sources for the GCRs, and the intrinsic com- 

bj plexity of this mechanism, the source spectral index and es- 

The measured Galactic cosmic-ray (OCR) fluxes at Earth re- pg^j^jjy jj, low-energy shape is not very well predicted (e.g., 

suit from a three-step journey: 1) the diff^usive shock acceleration Caprioli itd]l2010l) . Awaiting further progress along this line, 

(DS A) mechanism provides a source specti'um; 11) these paiti- ^„ j^^jj^g^t ^^ute to access the source spectrum is to start with the 

cles are then transported (by diff^usion, but also convection and top-of-atmosphere (TOA) fluxes and go back to the source spec- 

reacceleration) and also interact in the interstellar medium (ISM) ^.^ ^o do so, it is generally assumed that the steady state holds 

until they reach the Solar neighbourhood; 111) they enter the solar ^r the propagation, and that the first and second step ai-e inde- 

cavity where they decelerate due the eff^ect of solar modulation ^^^^^^^ xhe first assumption is known to fail at high enough en- 

(active for GCRs below a few tens of GeV/n). eigies, whereas the second one may only be approximate. In this 

The last step prevents direct measurements of low-energy study, given the success of simple steady-state diff^usion models 

(beyond a few hundreds of GeV/n) interstellar fluxes (IS). The foj. nuclear component, we follow the same route in order to 

first step can be investigated by means of semi-analytical or nu- draw some consti-aints on the source parameters. Note that this 

merical studies of the DS A mechanism. However, due to the va- paper does not consider primary electrons, whose specti-um may 
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depend on local sources abo ve a few tens of GeV (e.g. jBoularesI 
ll989HDelahave et al.ll2010l ). 

The propagation step is the focal point of many phenomeno- 
logical studies addressing the flux of secondary species (cre- 
ated by interactions of the primary species in the ISM and 
radiation fields of the Galaxy) including light nuclei, antipro- 
tons, positrons, radioactive isotopes and also gamma rays. The 
transport parameters are usually determined by fitting data on 
secondary-to-primary ratios of nuclei, as for example the B/C 
(boron-to-carbon) ratio. However, present data on B/C ratio, 
even when combined with CR radioactive isotope measure- 
ments, lead to a constrained b ut large range of allowed val- 
ues fo r these parameters (.Maurin et al.ll200ll 120101 : iPutze et al] 
I2OIOI) . The study of these secondary- to-primary ratios is — to 
first order — i nsensitive to the details of the source spectra (e.g., 
iMaurin et al.l |2002). Therefore, a common phenomenological 
approach is to first extract the transport parameters, then to fit 
the source spectra; however the source and transport parameters 
may be correlated (Putze et al. 2009). 

The importance of the primary fluxes and their ra tios was 
recognised a long time ago d Webber & Leznia^ll974l) . In this 
paper we reconsider their study, trying to answer the following 
questions: what phenomena shape the TOA and IS fluxes? Is it 
the source spectrum or the propagation step, and are there de- 
generacies between the two effects? To accuracy accuracy can 
we determine the low-energy source spectra, the spectral indices, 
and the source abundances? Are the source spectra universal or 
species-dependent? 

On the experimental side, accurate data are available up to 
few hundreds of GeV, whilst at higher energies data are less 
abundant and have large error bars and scatter between exper- 
iments. On the modelling side, we have a semi-analytical prop- 
agatio n model proved to work well with many GCR ob serv- 
ables (Mauri n et al.ll200Tl: lDonato"etani200ll I2001 120091) and 
an implementation of the Markov Chain Monte Carlo technique 
(MCMC) to derive the probabili ty density func tions of the anal- 
ysed model parameters (Putze et al.l2009ll2010l) . We take advan- 
tage of these to address the above questions. Our results are also 
supported by toy-model calculations, especially for the shape 
of ratio p/He. As abundant and accurate data are expected in 
the near future by the orbiting PAMELA experiment and forth- 
coming AMS-02 detector (to be installed on the International 
Space Station), such a study also aims at providing some guide- 
lines on how to tackle the information contained in the primary 
flux propagated spectra. We finally note that the recent ATIC- 
2 and CREAM-I measurements for proton and Helium data 
hint at a spectral change > TeV/n. This has consequences for 
the secondary production of y-ray s, antiprotons, and positrons 
dPonato & Serpicoll2010l : iLavalld 120101) . However, this occurs 
only for the high-energy part of thes e spectra. In particu lar, as 
our previous antiproton calculations (iDonato et alJl2001h rehed 
on a fi t to the data, the conclusions obtained in ^Don ato et al] 
(I2008h remain unchanged even with the new source spectra pro- 
vided in this study. 

The paper is organised as follows: Sect. |2] contains a short 
overview of the propagation scheme employed in the present 
study. In Sect. [3] we discuss the p and He data and whether they 
can provide any constraints on the transport parameters or if they 
can be fitted in any propagation configuration (e.g. with or with- 
out convection, with or without reacceleration). In Sect. |4l we 
seek for generic constraints on the source spectra (p. He, and C 
to Fe) comparing the values obtained in different configurations 
of propagation models. In Sect. |5] the origin for the observed 



shape for the ratio of primary species is outlined. Our conclu- 
sions and perspectives are given in Sect.|6l 



2. The propagation model 

The framework employed to calculate the fluxes is the dif- 
fu sion model with convectio n and reacceleration discussed 
in 'Maurin et al.' ( 20011 l2002h . updated and fully detailed in 
Putze et al. (2010). Here we only summarise the main features 
of the model. 

The Galaxy is shaped as a gaseous thin disk with half- 
thickness h = 0.1 kpc and an infi nite radial extension (ID model, 
as also used in Jones et al. 2001), hosting the interstellar medium 
and the stars, and surrounded by a thick halo for cosmic-ray 
transport whose half-height is L. 

Assuming a steady-state situation, the transport equation for 
the CR nucleus j can be written as 



(1) 



where the difl^erential density A^-' = N^(E, r) depends on the po- 
sition r in the Galaxy and on the energy (throughout the paper, 
E is the total energy, is the kinetic energy, Eic/„ the kinetic en- 
ergy per nucleon, and E/„ the total energy per nucleon). sums 
up the physics of the transport in the Galaxy, while S^' contains 
the source term. 



2. 1 . Transport parameters 

The operator X (we omit the superscript j) describes the diffu- 
sion K(r, E) and convection Vc{r) in the Galaxy, the decay rate 
^md{E) = l/(7To) for radioactive species, and the destruction 
rate Ti„f,]{r,E) - Y^is m ni5u{r)vcri„AE) on the interstellar mat- 
ter (ISM). It reads 



Zir, £) = -V • (KV) + V • Fc + r„d + r,„ei. 

The spatial diffusion coefficient is parametrised as 



(2) 



(3) 



where H - pcjZe is the rigidity of the particle, y6 is the velocity 
of the particle in units of c and tjj parameterizes the low-energy 
behaviour of diffusion. The nominal form is given by 777 = 1, 
which is just the inevitable effect of particle velocity on the dif- 
fusion rate. Hence this p term is always present, and other val- 
ues of rjj are relative to this. lPtuskin et al.l (i2006) argued that the 
form of the spatial diffusion coefficient could be modified at low 
energy, due to the possibility that the nonlinear MHD cascade 
sets the power-law spectrum of turbulence. Indeed. Mauri n et al] 
(2010) found that the value of this parameter was crucial for the 
determination of 6 given the current B/C data. The convective 
wind acts in the whole diffusive volume with a constant velocity 
Vc - ±Vcez pointing perpendicularly to the Galactic disk. The 
coefficients b and c account for the first and second order energy 
changes 



idE\ 

b(r,E) = ( — ) 

\ (j^f I ion, coul. 



m + Ek 



c{r,E) = P^xKp 
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Coulomb and ionisation losses add to possible energy gains 
due to reacceleration, described by the coefficient Kpp in mo- 
mentum space. The parameterisation for Kpp is taken from the 
model of minimal reacce l eration by the interstellar turbulence 
dOsbome & Ptuskinlll988l:ISeo & Ptuskiniri994i) : 



(4) 



where Va is the Alfvenic speed. 
2.2. Source parameters 

The source term S includes the initial spectrum at the sources 
and the secondary contributions (spallations of heavier n uclei). 
Acceleration models typically predict dQ/dp oc p " (e.g. JJonesI 
11994 ). which leads to dQ/dE oc p^" IP, where the low-energy 
behaviour is unknown. For future reference, we note that 



dQ _ dQ _\ dQ _\ dQ _ \ dQ 
Ie ~ dFk' a' dEk/„ ~ ^ ~dp ~ Zp"d^ 



(5) 



In this paper, we model the low-energy shape by adding one 
free parameter, tjs , active at low energy: 



dQ 

dE. 



k/n 



(6) 



where q is taken to be the normalisation for a differential energy 
per nucleon source spectrum. The reference low-energy shape 
corresponds to rjs - -1 (to have dQ/dp oc p^", i.e. a pure 
power- law). 

2.3. Free parameters of the model 

The present model contains of necessity several free parame- 
ters: the parameters in the transport sector {Kq, 6, Va, Vc, tjt}, 
the ones in the source term {q, rjs , a], and the halo size of the 
Galaxy L in the geometry sector We will see in the follow- 
ing that not all these parameters have the same relevance to the 
physics of primary cosmic nuclei, and we will therefore operate 
within a critical sub-range of parameters. In diffusion models, 
L cannot be solely determined from the B/C ratio because of 
the well-known degeneracy between Kq and L when only sta- 
ble species are considered. If not differently stated, we will work 
with the default values rjs - -I and 777 = 1, and the reference 
value L = 4 kpc. In most of the analyses we let free the source 
normalisation qi for each primary species. 

The low energy (< 10 GeV/n) charged particles are braked in 
the heliosphere by the solar wind, modulated with variations of 
the scale of the 11 -year cycle. We adopt the force-field approx- 
imation, which provides a simple analytical one-to-one corre- 
spondence between the modulated top-of-the atmosphere (TOA) 
and the demodulated interstellar (IS) fluxes, and whose only 
effective parameter is the modulation potential (p (GV). For a 
species j, the IS and TOA energies per nucleon are related by 
E^^ - EJ^^ -h O (O = Z/A X is the modulation parameter), and 
the fluxes by (p is the momentum) 



(7) 



The force-field method is an approximated solution to the 
diffusion equation for charged particles in the heliosphere. A 
more accurate treatment is far beyond the scope of our paper. 



which is centred on the source and diffusive processes active in 
much a wider energy range and responsible for strong effects of 
the primary spectra. 

In principle, the solar modulation potential could also be a 
free parameter of the study. Since the low-energy spectrum of 
the primary fluxes is determined mainly by the low-energy in- 
jection spectrum (see Eq. |6]l, the low-energy diffusion scheme 
(see Eq. Q), and the solar modulation (see Eq. |7]i, a non-trivial 
correlation is expected between these parameters. For instance, 
we find using the minuit minimisation routines (not showed) that 
there is a negative correlation between <p and rjs . This holds for 
all primary cosmic-ray fluxes studied in this work. Given that 
there is already a great number of degeneracies between the dif- 
ferent source parameters, which cannot be lifted with the current 
available data, we choose to fix the modulation potential to the 
values given by the exp eriments. A more d etailed study, in the 
spirit of the one done in iTrotta et al.l (1201 Ol) . will be undertaken 
in a future work. 



3. Analysis with free source and transport 
parameters 

In this Section, we use different data sets for proton and he- 
lium fluxe£] in order to determine which propagation models 
describe data and to try to set constraints on the free parame- 
ters of our model. The fitting procedure is based on the minuit 
routine, which minimises a x~ function. 

3.1. Data 

A first important consideration is the choice of the data used to 
constrain the models. The top panel in Fig. [T| shows the avail- 
able data for p and He fluxes. The abscissa is the kinetic energy 
per nucleon (Zst/n) and the ordinate i/r'^ x E^-^^, where i/r'^ is the 
IS, demodulated using the force-field approximation. The low- 
energy region (below 100 GeV/n) has been covered by many 
balloon-borne, Space-Shuttle based and satellite experiments, 
whereas above a few TeV/n the data come from several bal- 
loon long-exposure flights (accumulated over several flights in 
a decade). ATIC data cover the gap at a few TeV/n energy. The 
overall agreement between the data is fair, the scatter between 
the data being higher at high energy. 

For our analysis, the criterion is to select data samples cov- 
ering a broad energy range and consistent between experiments. 
We show a subset of demodulated data in Fig.|2]to illustrate the 
error bars and differences between the mos t recent and c onsis- 
tent sets of p and He data, name ly AMS-01 dAlcaraz et alj 



AMS Collaboration et all bOOOl) . BESS98 dS anuki et alJ 



200c ; 



200c ; 



Shikaze et al. 2007) and BESS-TeV (a.k.a. BESS02. lHaino et al.1 
2004; Sh ikaze et al.ll2007l) . For the proton flux, the AMS-01 and 
BESS98 data, both taken in 1998 in the same solar period, are 
consistent except at low energy. BESS-TeV data taken in 2002 
during a high level of solar activity show a different behaviour at 
low and intermediate energies. Note that usually the solar mod- 
ulation level is obtained by fittin g and a simple tw o-parameter 
proton spectrum to the data (S hikaze et alj|2007l) . This is the 
standard lore in the field, although it is expected to give a biased 



' The isotopic separation is not always achieved, so that many data 
actually correspond to 'H+^H for the proton flux and ^He+''He for 
the helium flux, where and 'He come from secondary contributions 
only. This amounts to a < 10% error at low energy, which is contained 
in the error bars. 
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10^ — 




SOKOL on COSMOS (1965-1966) 
MUBEE (1975-19S7) 
JACEE (1979-1995) 
Ichimura et al. (1987-1991) 
IMAX(1992) 
BESS (1993) 
CAPRICE (1994) 
RUNJOB (1995-1999) 
RICH-II (1997) 



BESS (1997) 
AMS (1998} 
BESS (1998) 
CAPRICE (1998) 
BESS (1999) 
BESS (2000) 
BESS (2002) 
ATIC-2 (2002) 
CREAM-I (2004) 



the series of BESS ballo on flig hts BESS9 3 dWang et al.l|200 ; 
BESS97 (Shikaze et al] l2007l) . BESS98 dSan uki et al. 2000'- 

200 7|). BESSOO 



IShikaze et al. 2007.) BESS99 dShikaze etal 
(Shikaze etai rl2007[). and BESS02 (BESS-TeV, 
[200 4: Shikaze et al.l l2007h . for which = 700 MV, 491 MV, 
591 MV, 658 MV, 1300 MV, 1109 MV res pectively. The inter - 



Haino et al] 



mediate an d high- energy data are: ATI C-2 dPanov et al. 



oil), JACEE (f A sakimori e t al 
RICH-II 



200' 



CREAM-I dAhn et alJl20 

MUBEE (IZatsepin et al.l 

RUNJOB JOerbina et al.' 2005), SOKOL divanenko et alJIT993l 
and Jchimura et al.. (.1993.) 



1993), 



1998! 

(Diehl et al. 2003: 




K AMSOl (*=600 MV) 
o BESS-TeV (<1>=1 109 MV) 
o BESS98 (0=591 MV) 

I 1 ^ 



1 10 10" 10 10' 10' 10° 

Ek/n [GeV/n] 

Fig. 1. Demodulated p and He data {xE^I^^) as a func- 
tion of Ekin- The low-energy data and the associated mod- 9^ 
ulation param eter are: AMS- 01 (Alcaraz et al. 2000 for H '■^ 
and IAMS CoUaboratiori et all F2000| for He) with = S 
600 MV; CA PRICE94 d Boezioet all [19991) and CAPRICE98 \ 
dBoezioetal.ll2003h with 6 ^ 65 MV and 600 MV respec- 
tively; IMAX ("MennetaLl 12000^ with ^ 750 MV ; and 



Ek/n [GeV/n] 




AMSOl {<1>=300MV) 
BESS-TeV (*=55-'i MV) 
BESS98 {■J>=295 MV) 



10^ 

Ek/n [GeV/n] 



Fig. 2. Demodulated p (top panel) and He (bottom panel) flux 
xEfl^ as a function of Ekjn for AMS-01 (black stars), BESS98 
(open red squares) and BESS-TeV (open blue circles). 



modulation level (for example, we do not know the true interstel- 
lar proton spectrum, there is the problem of polarity in the solar 
magnetic field, etc). The use of demodulation has the same phys- 
ical basis as modulating in the same scheme (the force-field in 
our case, see Eq.|7]i. In Figs. 1 and 2 we have chosen to demodu- 
late data to compare all existing data on the same foot. Since the 
solar modulation is effective up to few GeV/n, only the very low 
side of the figures are affected by the demodulation procedure. 
The goal of this paper is not to deal with these issues, but sim- 
ply the fact that the proton data are already inconsistent among 
themselves implies that we may expect inconsistencies in the fit- 
ted models. 



3.2. Pure diffusive transport 

The first step is to test a minimal model containing only acceler- 
ation and plain diffusion, as well as nuclear reactions and elec- 
tromagnetic energy losses, but without convection and reaccel- 
eration (y„ = = 0). 



Main degeneracies In the pure diffusion model, the flux of any 
primary species at high energy can be approximated b}Q 



Q{E) q 
w(E) oc oc — ■ E 



(a+S) 



(8) 



This formula shows two degeneracies between the source and 
transport parameters: the first one is in the normalisation q/Ko, 
the second one is in the total spectral index a + 6. 

We start with a minimisation procedure setting the free pa- 
rameters 6, a and He (in order to break the degeneracy between 
q and Kq the latter is set to 0.0048 kpc^ Myr"'). The results for 
protons and helium are presented in two different columns in 
Table [T] The Xmm/d o f values are very small for a number of 
cases, indicating a possible over-fitting of the data. The values 
of the best-fit parameters are not reported since they are not rel- 
evant at this stage of the analysis. For the different sets of data, 
the values of both a and 6 vary from almost any value between 
and 2.8 (not shown in the Table), but the sum of them is close to 
3.0 for p and 2.8 for He. The first three lines show that a fit to p 
and He on the AMS or BESS data is always possible in a simple 
diffusion scheme (although it provides unphysical values for a 



- Throughout the paper, the quantity if/ denotes the differential flux 
in kinetic energy per nucleon, i.e. dif//dEk/„. The notation dif//d'R used 
in Sect. |5] is the only place where it will refer to the differential flux in 
rigidity. 
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Table 1. Best-fit to p and He data for pure diffusive transport 
(experiment, number of data points, best ;\f^/d.o.f., and best-fit 
a + 6 value). 



Data 




P 






He 






# 


d.o.f. 


(y '\~ S 


# 


d.o.f. 


CK + (5 


AMS-01 (1) 


28 


0.19 


2.99 


31 


0.72 


2.81 


BESS98 (2) 


41 


0.39 


2.99 


36 


0.36 


2.77 


BESS-TeV (3) 


47 


0.72 


2.89 


40 


0.80 


2.80 


(1+2) 


69 


0.66 


3.00 


67 


1.75 


2.81 


(1 -H 2-^3) 


116 


1.95 


2.99 


107 


2.52 


2.81 


All LEt 


304 


6.22 


2.96 


287 


2.77 


2.79 


ATIC-2 


15 


22.46 


2.75 


15 


37.4 


2.81 


All HE* 


59 


2.58 


2.87 


39 


4.18 


2.85 



IMAX92, BESS93, CAPRICE94, BESS97, AMS-01, BESS98, CAPRICE98, BESS99, 
BESSOO, BESS02. 
* CREAM, SOKOL, MUBEE, JACEE, Ichimuia et al., RUNJOB. 



and (5). The fourth line shows the co mbined analysis of A MS- 
01 (lAlcaraz et al.ll2000h and BESS98 dSanuki et al.ll2000h data. 
They were collected in the same year (1998) — which may help 
reduce the systematics due to solar wind modelling — and span 



nearly the same energy range. The x' 



min/d.o.f. 



1.75 for the fit of 



combined He data, compared to the corresponding ^ ^ val- 
ues of 0.72 and 0.36 for the separate data shows inconsistencies 
among the data sets, as underlined in Sect. 13.11 When combin- 
ing the three experiments (fifth line), this is even more visible, 
also for protons (the best-fit slope is unaffected). We have fitted 
also all the available data in the low energy (sixth line) and in 
the highest energy range (last line). The agreement among the 
data sets is poor, both for protons and helium, as akeady visi- 
ble in Fig.[T] ATIC data, which connect the low and high energy 
sectors, are badly fit by any pure diffusive model. 

We may naively interpret the results of Table [T] as showing 
that p and He data can be well accommodated in any purely dif- 
fusive transport models, in the low energy range (< 100 GeV/n), 
for each experimental data set taken separately. But such models 
lead to unphysical values for 6 and a. So it could also mean that 
the hypothesis of a standard source spectrum (i.e., dQ/dp oc " 
when setting rjs - -1) and a standard propagation scheme 
(rjT = 1) is unsupported by the data, or that additional effects 
(e.g., convection and/or reacceleration) are required to match the 
data. Before resolving this issue, we go further with the compar- 
ison of the approximate formulae and the full calculation. 



Inelastic interaction: a link between a, 6 and Kq. Explicitly, 
the effect of the catastrophic losses at low energy, Eq. ([8]l gives, 
for a ID model. 



HE) = T- 



hL~ + "ismo-v 



(9) 



For 6 fixed, Eq. Q implies a coiTelation between Kq, a and 6, 
given some primary data. Indeed, as Kq decreases, the inelastic 
interaction term nisMCt^' becomes more efficient in the denom- 
inator of Eq. (|9]l. The effect of species destruction is therefore 
more pronounced at low energy. Fixing 6 and going to small Kq 
we expect that, in order to balance the increased destruction rate, 
the numerator compensates but increasing a. If Kq is fixed and 
small, so that inelastic interactions can dominate, the same flux 
of protons (or helium) can be obtained with a larger a + 6. 

This effect is confirmed by the numerical results, as seen in 
Fig. [3] For each point in the K() - 5 plane, we plot a + 5 for the 




Fig. 3. Surfaces of a + y for the best-fit models in the plane 
Kq - 6 (free parameters are a and ^p) on AMS-01 proton data. 
The colour code (from light to darker shades) for the contours 
superimposed on top of each graph correspond to (a + y) - {2.8, 
2.9,3.0,3.1,3.2}. 

best-fit model on AMS-01 proton data (the free parameters are 
a and q^). We checked (not shown) that similar values for a + 6 
are obtained when the fits is performed on other proton fluxes 
(BESS98 or BESS-TeV), or for other species (He or a combined 
fit p-i-He). We can see from the figure that, for any fixed 6 value, 
the data require higher a + S while Kq decreases, due to the in- 
creasing importance of the destruction rate. This effect is less 
pronounced for small values of the diffusion coefficient slope, 
namely when S is close to 0.2-0.3. 

Asymptotic behaviour (or why a + 6 > ydata) Even for Hght 
primary species such as protons, which suffer less from destruc- 
tion in the ISM, the asymptotic purely diffusive regime is not 
reached. If we fit a primary flux with (as is usually done in the 
literature) 

il/{E) oc E'^"'" , 

then we are bound to have 

a + 6 = rasympt ^ Tdata ■ 

This implies that caution is in order whenever we wish to com- 
pare the result of studie s fitting the pr opagated fluxes with a 
power-law function (e.g.. lShikaze et al. 2007) to those (such as 
this one) fitting directly the source spectrum. Catastrophic, but 
also continuous losses flatten the propagated spectrum below < 
few tens of GeV/n energies. The inequality yasympt <: Tdata is also 
valid for conve ction and/or reacceleration, as derive d from an 
analysis of B/C (iPutze et al.ll2"oToHMaurin et al.ll2010l) . 



Simultaneous fit of p and He to lift the a + S degener- 
acy? As the residence time — hence the destruction rate of any 
species — depends on the energy through the transport parame- 
ter 6 (and not ona + 6), the different inelastic cross-sections for 
each species (cr'' , ~ 30 mb and cr^'^, ~ 90 mb) leave different 

^ inel inel ^ 

imprints on the corresponding p and He low-energy spectra. This 
is expected, to some degree, to lift the degeneracy ona + S when 
using a combined fit to various primary species. 

Figure |4] shows the A'min/d o f contours in the Kq - 6 plane 
for the separate fits of p (first row). He (second row) and for the 
combined fit p-i-He (last row), for the different sets of data used 
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Fig. 4. Surfaces of x^-^^^/^^^f (best-fit models) in the plane Kq - 6. First row: fit performed on p data (free parameters a and ^p). 
Second row fit performed on He data (free parameters a and ^He)- Third row: fit performed simultaneously on p and He data (free 
parameters a, q-p and ^He)- The columns from left to right correspond respectively to AMS-01, BESS88 and BESS-TeV data. The 



colour code for the contour plots {^^o-x]^h^ia of) superimposed on each graph is (from white to darker shades) x\ 
5, 10). 



min/d.o.f. 



{0.5, 1, 2, 



before. The ranges chosen for Kq and 6 correspond to extreme 
but not impossible values of these parame ters that can accom - 
modate the secondary-to-primary B/C ratio (iMaurin et al I20T0I) . 
The top-left plot shows the strong degeneracy of a and 6 (for 
AMS-01 data), as almost any configuration is acceptable. The 
top-right plot shows that for other data (BESS-TeV) no good fit 
can be achieved ixlai^ii of > 2) in the selected Kq-6 region: the 
good fits occur for unrealistic 6 only. We note here that there is 
no inconsistency with the results shown in Table [T] which have 
been obtained by spanning larger ranges for the free parameters 
including unphysical values. The second row of Fig.|4]shows the 
fits to He data. All the experiments tend to prefer large Kq and 
large 5 (both unrealistic if we demand ^ ^ ~ 1), but the;\f^ 
does not change significantly, meaning that no particular class of 
models is selected by the helium data. 

The third row of Fig. |4]is the combined p-i-He fit. The re- 
sulting ^ ^ surfaces corresponds to a trade-off between the 
best-fit for p and He. The best-fit 6 (not shown in the figure) still 
falls in region of 5 > 1, and the degeneracy a + 6\s, not lifted as 
no specific value for 5 is preferred. 



The role of rfs and rjj- We have introduced the possibility to 
have a non-standard low-energy diffusion coefficient by means 
of the parameter 777, as given by Eq. The low-energy shape 



of the source spectrum is driven by the parameter rjs , see Eq. (|6]l. 
If for illustration we neglect the nuclear interactions, we obtain 



QiE) 
K{E) 



(10) 



where the extra p factor comes from the vjAn in front of Eq. 

The parameters rjs and tjj introduce a similar shape cor- 
rection at the lowest energies. If it were not for energy losses 
and inelastic reactions (see next section), only the quantity 
rjs - TjT would be expected to be constrained. We fit AMS-01 
and BESS98 data, as well as all the proton data, with a, rjs, S, 
and as free parameters. The best;^'^ is slightly smaller than the 
one obtained with only a, 6, and free (and tjs - -1). Similar 
results are achieved when the acceleration scheme is fixed to the 
standard lore (775 = -1) and the fourth free parameter is 777. 
However, the corresponding values for 6 and a are again unsup- 
ported. 

A further degeneracy can be produced by solar modulation, 
whose action is a decrease of the flux with the increase of the 
solar wind strength (parameter (p). On the other hand, the TOA 
flux increases with (p when rjs < tjj (in the standard scenario 
rjs - -1 and 777- = 1), so that some compensation of the solar 
modulation can be produced. 
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3.3. Summary for the transport parameter constraints from p 
and He data 

The existing proton and hehum data are unable to select any par- 
ticular propagation model. This is consistent with the fact that 
the transport and source parameters are degenerate, as shown 
from simple arguments in our toy formulae. The best present 
data (AMS-01, BESS98 and BESS-TeV) can be quite well repro- 
duced by solar modulated pure diffusive transport for instance, 
but they favour unphysical values of the source and transport 
parameters. Hence other ingredients are required. This could be 
a modification of the low-energy source spectrum or diffusion 
coefficient, or the addition of convection and/or reacceleration 
(which can also accommodate the current data), or an improve- 
ment of the calculation of the solar modulation effect. However, 
given the present data and the physics of primary spectra, such 
an approach is bound to fail: the increase of the parameter space 
merely brings new degeneracies. Moreover, most of the param- 
eter space is already ruled out by the constraints from the B/C 
ratio. To go further, we thus have to restrict the parameter space 
to the source parameter space only (and use some prior on the 
propagation parameters). 

4. Analysis with fixed transport parameters 

In the previous sections we have shown that present data on 
primary fluxes alone cannot constrain significantly the trans- 
port parameters. The next natural strategy is to fit simultane- 
ously primary and secondary species. However, as emphasised 
in Putze et al.] (l20()9l) . the large body of data for primary species 
drives the fit away from the best-fit regions of the B/C ratio. 
Actually, the standard lore is to fix the transport parameters to 
their best-fit value, and then constrain the source parameters. But 
the latter values are then biasecfl. We nevertheless follow this ap- 
proach, but we repeat the analysis on several possible transport 
configurations. This allows us to derive explicitly the systematic 
effects of the source parameters arising from this bias. 

In this section, we first gather several sets of transport param- 
eters shown to be consistent with B/C data (Sect. l4Tb . We then 
fit the source parameters for p and He, the best-measured pri- 
mary fluxes to date (Sect. l4!2b . We repeat the analysis for other 
primary species, to inspect the universality of the source slopes 
and obtain their relative source abundances (Sect. l431 l. 

4. 1 . Transport parameters consistent witli B/C 

The transport parameters are usually constrained from 
secondary-to-primary ratios (e.g. B/C). In the literature, 
various classes of models have been used, leading to very 
different values of their respective best-fit parameters (see, e.g., 
IStrong et al.l |2007| for a review). For instance, a model with 
diffusion + reacceleration is ch aracterised by a best- fit propaga- 
tion slope 6 ~ 0.3 - 0.4 (e.g., iLionetto et d]l2005l). leakv-box 
inspir ed models points to 6 ^ 0.5 - 0.6 (e.g., IWebber et al.l 
120031; IPutze et al.. ,2009.) . whereas diffusion + convection 
models (w/wo reacceleration) points to 6 ~ 0.75 - 0.85 (e.g., 
iMaurin etal ] |200Tl IPutze etal]l20Tol) . 

This sensitivity to the CR transport mode (pure diffu- 
sion, w/wo convection, w/wo reacceleration) is discussed in 

^ IPutze et al. (2009) show that the values of the source-spectrum pa- 
rameters (slope and abundances) are positively correlated among them- 
selves and with the reacceleration strength, but are negatively correlated 
with the other propagation parameters. 



Table 2. Best-fit transport parameters for various configurations 
of the diffusion model (fitted on B/C data). 



Model 


rjT 


K'^''' X 10- 
(kpc^ Myr" 


) 


ybest 

(km's-') 


l/besl 

(km s" 


^7d.o.f 

') 


II 


1. 


9.76 


0.23 




73.2 


4.73 


III 


1. 


0.48 


0.86 


18.8 


38.0 


1.47 


I/O* 


-2.6 


2.05 


0.61 


0. 




3.29 


III/II 


-1.3 


3.16 


0.51 


0. 


45.4 


2.26 



Note 1. These results were obtained for input ingredients described in 
IMaurin et al] <2010l) . and for L = 4 kpc. 

iMaurin et al ] (120101) . Although the best-fit model is one with 
both convection and reacceleration, it predicts 6 ~ 0.8, a value 
quite high compared with theoretical expectations (1/3 for a 
Kol mogorov spectrum of turbule nce, and 1/2 for Kra ichnain, 
e.g.. lStrong et al.l l2007^. Following Maurin et al] (l2010l) . we use 
below four configurations of the diffusion model covering a large 
but plausible range for the transport parameters. These models 
along with their best-fit parameters are reproduced in Tab.|2] 

- Model II is with reacceleration only; 

- Model III is with convection and reacceleration; 

- Model I/O is with a low-energy upturn of the diffusion coef- 
ficient (riT < 0); 

- Model III/II is as I/O, but with reacceleration. 

The first two models correspond to the best-fit parameters for a 
standard spatial diffusion coefficient [i.e. rjr set to 1, see Eq. ^]. 
The last two lines correspond to a modified diffusion scheme: 
negative values of rjr are associated to an u pturn of the diffusion 
coefficient at low energy (iPtuskin et al .'2006). These two models 
are respectively termed I/O and III/II because both allow some 
convection, but both favour Vj"^'*' - 0. As shown in Fig. 7 of 
IMaurin et al.l (l2010l) . these models fit reasonably well the B/C 
data. 

4.2. Constraints on p and He source parameters 
4.2.1. Generalities 

On the one hand, the low-energy shape of the source spectra 
are not well known theoretically. They result from the diffu- 
sive shock a cceleration mechanisms at pla y in supernova (e.g . 
lDrurv"l983') or super-bubble shocks (e.g. Iperrand et al.ll20oi 
Ferrand & Marcowith 2010). Power-laws close to -2 (in energy 
space) are predicted at high energy, but there is still no agre e- 
ment about the low-energy spectrum (e.g. CaprioU et alj|2010|). 

On the other hand, we have access only to propagated spec- 
tra, where effects such as destruction on the ISM, energy losses. 
Galactic winds and reacceleration change the energy spectra up 
to a few tens of GeV/n. This is the route followed in this sec- 
tion, where we try to constrain the source parameters from a fit 
to the propagated fluxes. However, the shape of the flux at low- 
energy is only an extrapolation since the low-energy interstellar 
spectrum is masked by solar modulation effects. 

Note that the low-energy (below 100 MeV) IS spectrum 
can be indirectly constrained due to its interaction with the in- 
terstellar medium. In this appr oach, the IS flux is based on 
empirical fits to the data (e.g., Herbst et"ani2010l) . and its ex- 
trapolation at low energy is used to calculate, e.g., the ionisa- 
tion o f the ISM (feebber 1987; Nath & Biermann 1994; Webb^ 
Il998h or and of molecular clouds (.Padovani et al.n2009.) . or the 
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Fig. 5. Left panels: PDF of the source slope a for p (unfilled 
histograms) and He (hatched histograms). Right panels: PDF for 
^He - fp- The colour code corresponds to the three experimental 
data used: AMS-01 (solid black line), BESS98 (dashed red lines) 
and BESS-TeV (dash-dotted lines). 



LiBeB Galactic enrichm e nt and production jGilmore et al.l 19921: 
iNath & Biermannl 119941: lLemoineet"al]|1998h . Some of these 
studies find an increase in the low-energy spectrum, others a de- 
crease (with respect to a pure power law). Actually, data from 
the Voyager 1 & 2 spacecraft near the heliospheric termination 
shock could also be helpful for such studies, as they are close 
to interstellar cond itions, their level of modulation being »; 60 
MV (IWebber et al.ll2008L IWebber & HigbiSIIOOl . However, it 
has been argued recently that anomalous cosmic rays could con- 
tribute to an important fraction of the proton spectrum below 
300 MeV (Schereret al. 2008). For this reason, we do not in- 
clude Voyager data in our fits, and will only compare them to 
the best-fit spectra (based on the other data) at the end of this 
section. 



4.2.2. Results 

Due to the lack of robust information about the low-energy spec- 
trum, we choose to rely on a simple parametrisation allowing 
for an increase or decrease at low-energy, as given by Eq. (|6]l, 
i.e. Qe^JE) - q ■ fi^'^ ■ . For each configuration given in 
Table|3[— i.e., for a given choice of the transport parameters /To, 
5, rjT, Va, and Vc — we then use the MCMC technique to get the 
probability density function (PDFfl of the three source param- 
eters qi, Tj'g and a, (where / is either p or He). The three data 
sets on which we base the analysis are AMS-01, BESS98 and 
BESS-TeV (see Sect.©. 

PDF of Q-p, ffHe, and an^-a^. The two left panels of Fig.|5]show 
the PDF of ffp (hatched histograms) and ane (empty histograms). 



The MCMC techniq ue and the in t erface with the propagation code 
is expla ined in deta i ls in [ Putze et al.l ( l2009h . A shorter introduction is 
given in lPutze et all 
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Fig. 6. Most-likely value (symbols), 68% and 95% CIs for p 
(filled symbols) and He (open symbols) for the four propagation 
configurations gathered in Tab. |2] Top panel: spectral index a. 
Bottom panel: low-energy source parameter rfs [see Eq.(|6]l]. The 
grey line tjs - - \ coiTesponds to the value for which the source 
spectrum is a pure power-law in rigidity (i.e. dQ/dfi oc 'Rr"). 



whereas the two right panels show the PDF of ane - Q'p to visu- 
ally inspect any discrepant spectral index for the two species. 
For AMS-data (solid black lines), both Model II (reaccelera- 
tion, 6 - 0.23) and Model III (reacceleration and convection, 
6 - 0.86) show a very good agreement between their p and 
He spectral index, with respectively au 2.45 and am ~ 2.3. 
There are significant differences for BESS98 (red dashed lines) 
and BESS-TeV (blue dash-dotted lines) data: first, the match be- 
tween (THe and Op is not as good as for AMS-01, yet aue - re- 
mains marginally consistent with 0. The plots on the right panels 
and the width of the PDF tell us that the data current precision 
does not allow us to separate differences < 0.1 in the spectral 
indices. 

The 68% and 95% CIs on the spectral indices for the four 
transport configurations of Table |2] and the three sets of data are 
shown in the top panel of Fig.|6] From a quick visual inspection, 
the following trends are found: 

- for any given data set and species (p or He), the spread in 
the source slopes is Up - ofee ^ 0.2, regardless of the model 
(values in the same column in Table|3]i; 

- for any given model, the typical spread in a when fitting dif- 
ferent data sets (values in the same row in Table|3ll is ~ 0.05 
for ap and a; 0.1 for aue- This is larger than the eiTors ex- 
tracted from the minuit minimisation routine, which gives a 
statistical uncertainty Op - aee S 0.01 (not shown); 
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- the spectral indices obtained from BESS-TeV data are sys- 
tematically larger and marginally incompatible with those 
found for AMS-01 and BESS98. This may be related to the 
systematically higher value obtained for rjs (see below). 

- model II (reacceleration only) gives larger spectral indices, 
inconsistent with the values found for the three other models. 
This is not unexpected as it has the smallest S of all models 
(considered in Table |2]i. 

A scatter of ~ 0.2 is thus attributed to the fact that we do not 
know which model is best, and a scatter ~ 0.1 because of sys- 
tematics in the data. 

Low energy and confidence intervals (CIs) on rjs The 68% 
and 95% CIs on the parameter ris controlling the low-energy be- 
haviour of the source spectrunfl are shown for the same mod- 
els/data in the bottom panel of Fig. |6] BESS-TeV data be- 
ing at slightly higher energy than AMS-01 and BESS98, its 
source spectrum low-energy parameter rjs is less constrained. 
Otherwise, the p and He t/s point to fairly similar values for any 
given propagation configuration. However, this value depends 
on the model chosen: the reacceleration model (II) and convec- 
tion/reacceleration model (III) both favour ris ~ 1, whereas rjs 
is close to -2 for Model I/O and -1.5 for Model III/II. The latter 
value is consistent with a source spectrum being a pure power- 
law in rigidity, whereas the former value implies a flattening at 
low energy. This is understood if we inspect the quantity 775-/77-, 
appearing in Eq. (fTOl i: for models II and III that have jjt - I, this 
give Tjs -r]T ~ 0. For models I/O and III/II that have respectively 
rjT = -2.6 and -1.3, this gives 775 -777- ^ 0.6 and -0.2. So it seems 
that the constraint Tjs-rjT ~ should be met for any propagation 
model. 

Source abundances <7, The scatter is quite large when all the 
different models/data are considered. The absolute values are not 
meaningful since they depends on the choice of L that is arbi- 
trary set to 4 kpc in this analysis. We nevertheless note that the 
ratio qHs/lp falls in the range 0.3-0.6 (not shown), with a typical 
spread of « 0. 1 - 0.2 for the PDF 

Spectra, goodness of fit, and liigli-energy asymptotic regime 
The data along with the best-fit spectra for all models are shown 
in Fig. I2] An eye inspection shows a good match to the data. 
More precisely, the ^ values given in Table [3]tell us that 
the fit to the data is very good for BESS98 and BESS-TeV, 
but not satisfactory for most of the models with AMS-01 data 
(for which the first and last two bins are not well reproduced 
given their small error bars). We remark that the spread in 6 is 
larger than the spread in a (see above). Hence, the smaller 6, 
the smaller 'Xasympt(= o: + 6). This is consistent with the same or- 
dering for all species of the propagated spectra — from larger to 
smaller yasympt — seen on Fig. Q The top curve is always Model 
II {6 = 0.23, solid Hnes), going down to Model III/II (6 = 0.51, 
long dash-dotted lines). Model I/O {S - 0.61, dashed lines), and 
then Model III {6 - 0.86, short dash-dotted lines) for the bottom 
curve. This emphasises that more accurate data in the high en- 



^ We underline that for IS fluxes, the lowest energy data points are at 
~ 0.8 GeV/n for p and ~ 0.4 GeV/n for He, see Fig.|2] This corresponds 
to /3p ~ 0.8 (fine ~ 0.7), so that having a pre-factor fH^ amounts to a 
difference of ~ 1/2 (~ 1/3 for He) with respect to the case rjs = for 
this low energy point. 




kinetic energy [GeV/n] 




Helium X 10" 



liinetic energy [GeV/n] 




Fig. 7. TOA (modulated) fluxes (times E^-^^) as a function of the 
kinetic energy per nucleon, for p and He. The symbols are black 
circles for AMS-01, red squares for BESS-98, and blue triangles 
for BESS-TeV. The curves correspond (for the four models of 
Table|2]i to the best-fit spectra obtained by a fit on the AMS-01 
data (top panel), BESS98 data (middle panel) and BESS-TeV 
data (bottom panel). 

ergy regime (Te V-PeV) are needed to better constrain the asymp- 
totic behaviour 

Envelopes on IS fluxes and consistency with low-energy 
Voyager data Finally, the three panels of Fig. [8] shows the en- 
velopes on the IS fluxes (obtained from the 95% CIs on the 
parameters). The flux is extrapolated down to an IS energy of 
~ 0.1 GeV/n, where th e demodulated Voyager energies fall 
(IWebber & Higbid l2009fl On the same plots are shown the 
demodulated AMS-01, BESS98, BESS-TeV, and Voyager data. 
Showing here the demodulated envelopes underlines the spread 

^ The estimated modulation parameter is 60 MV for these data 
iWebber & Higbiell2b09i) . 
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Fig. 8. 95% CL envelopes for the proton and helium fluxes for the four propagation configurations of models of Table |2l The three 
panels correspond respectively to the result of the MCMC analysis on AMS-01 data (left panel), BESS98 data (middle panel), 
and BESS-TeV (right panel). For the sake of comparison, the IS (demodulated) AMS-01 (black circles), BESS98 (red squares), 
BESS-TeV (blue triangles), along with the Voyager data (stars). 



of the different models at low energies and the possible con- 
straints which can be obtained by low-energy data, such as from 
the Voyager spacecrafts. The low-energy envelopes obtained 
with AMS-01 and BESS98 data at Voyager energies are hardly 
consistent with each other. BESS-TeV envelopes are more sat- 
isfactory in that respect. The Voyager data for helium is well 
reproduced. This is possibly related to the use of the force-field 
approxi mation which is known to fail for the very low-energy 
protons (|Perkd[r987h . For the envelopes based on BESS-TeV 
data (right panel), almost all models are allowed, except perhaps 
model III (standard diffusion with convection and reaccelera- 
tion). On the other hand, from the envelopes from AMS-01 and 
BESS98 data, the modified diffusion scheme models (I/O and 
III/II) are disfavoured by Voyager data (that were not included 
in the fit). Again, it is difficult to conclude given the inconsisten- 
cies between the various data sets, but such plots clearly show 
the potential of future analysis to access the low-energy source 
spectrum (using Voyager data closer to the IS state and/or more 
accurate low-energy data from PAMELA and AMS-02). 



4.3. Constraints on heavier primary species 
4.3.1. Preamble 

Heavier primary species (from C to Fe) are less abundant and 
thus more difficult to measure than p and He. As a result, the 
spread in their measurements and their error bars are larger than 
those for the p and He fluxes. But they still provide some use- 
ful information on cosmic -ray propagation and sources. Indeed, 
the heavier the species, the larger its destructive rate. Hence, the 
universality of the source spectrum can be checked against the 
above effect (which is species dependent). 

In this section, we repeat the analysis performed on p and He 
for the C, O, Ne, Mg, Si, S, Ar, Ca and Fe elements. These ele- 
ments are almost all completely dominated by the primary con- 
tribution, except for S and Ar that receive a ~ 20% secondary 
contribution. To speed up the calculation, but still take into ac- 
count this contribution, we separate the nuclei to propagate in 
three families: '^C-^"Si, ^^S-'^^Ca, and ^'*Fe-*'^Ni. For instance, 
the nitrogen — a mixture in comparable amount of primary and 
secondary contribution — is properly calculated in this approach. 
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Table 3. Best-fit spectral index a for p and He fit and associ- 



ated Xmin/d o f (models/data coiTespond to those shown in Fig. |5] 
and the number of data for the fit is roughly the same for each 
species). 



Model/Data 


AMS-01 


BESS98 


BESS-TeV 




' a.o.f 













- Protons — 



n 


2.46 


2.14 


2.45 


0.70 


2.54 


0.48 


m 


2.30 


3.72 


2.30 


1.75 


2.36 


2.07 


111/11 


2.24 


2.05 


2.23 


0.73 


2.30 


0.93 


I/O 


2.28 


0.14 


2.28 


0.18 


2.32 


0.56 






— Helium — 








n 


2.45 


3.54 


2.41 


0.69 


2.53 


0.34 


m 


2.30 


2.85 


2.32 


0.41 


2.39 


0.40 


111/11 


2.27 


2.02 


2.26 


0.33 


2.36 


0.25 


I/O 


2.27 


1.00 


2.30 


0.21 


2.36 


0.24 



and could be in principle used in this study. However, we prefer 
to focus on the pure primary contribution to simplify the analysis 
and the discussion (a full analysis of all nuclei is left to a future 
study). 

Such a study complements a nd extends the analysis per- 
formed by the HEAP- 3 group (Eneelmann et al] Il990l) . the 
Ulysses group ("Duver nois & Thavei 19961) . and the TRACER 
group (Ave et al. 2009), in which only one propagation model, 
a universal source spectral index a for all species, and a sin- 
gle experiment was considered. Below, ff,, qi and rj'^ are free 
parameters for each primary species, which allows us to i) test 
the universality of the source spectra, ii) take into account the 
correla tions between the normalisation qi and the spectral in- 
dex a, dPutze et alj|2009l) . and iii) inspect the systematic spread 
on the source parameters as several configurations of the diffu- 
sion model are taken. The goal is to get more robust results (as 
more potential sources of uncertainties are taken into account). 
Besides, the MCMC technique is again helpful in providing a 
sound statistical estimate of the error bars on the source param- 
eters. 

We restrict ou r analysis to the HEAO-3 (lEngelrnann et al 



1990), TRACER (lAve et alJ l2008l) and CREAM-II (lAhn et al 
2009) data, to keep only the data covering as large as possible 
an energy region, and also to avoid very low-energy data which 
are more sensitive to solar modulation. 



4.3.2. Results 

We first show the fits to the data in Figs.|9l[T0l andfTTI They all 
show the same sets of data, but each figure corresponds to a fit 
to a single experiment (respectively, HEAO-3, CREAM-II and 
TRACER). Note that each species has three free parameters a 
(slope) , q (normalisation), and rjs (low-energy behaviour). 

Fit to HEAO-3 data In Fig. |9] (i.e. fit to HEAO-3), the four 
propagation configurations of Table |2] lead to the same shape 
at low energy. This is not surprising since this is where the bulk 
of HEAO-3 data lies. The high-energy asymptotic behaviour is 
then influenced by the value of the diffusion slope 6, as is the 
case for p and He. Similarly, the spread in 6 is larger than the 
spread in a (see below). It means that the smaller the value of 
6, the smaller yasympt - a + 6, so the same ordering according 




10' lO" 10' 

kinetic energy [GeV/n] 



Fig. 9. TOA (modulated) fluxes (times Ef^'^^) as a function of 
the kinetic energy per nucleon, for primary s pecies from C to 
Fe. T he symbols are black circles for HEAO-3 (lEngelmann et"an 



119901) . red squares for TRACER (Av e et al.11200 8') and blue tri- 
angles for CREAM-II ( Ahn et al. 2009|). The curves correspond 
(for the four models of Table |2| to the best-fit spectra obtained 
by a fit on the HEAO-3 data only. 



to the yasympt valuc of the model is seen at high energy: Model 
II on top (largest yasympt), then Model III and I/O, and Model 
III/II at bottom. An eye inspection shows that Model I/O and 
Model III/II are the ones in best agreement with the higher en- 
ergy data (CREAM-II and TRACER). 

Fit to CREAIVI-II and TRACER data Fig. [TOl shows the re- 
sulting best-fit spectra for the same models, but now fitted to 
CREAM-II data only. The data being at higher energy, the pa- 
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Fig. 10. Same as in Fig.|9] but the source spectra are now fitted 
on the CREAM-II data only (no published data for S, Ar, and 
Ca). 
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Fig. 11. Same as in Fig.|9] but the source spectra are now fitted 
on the TRACER data only. 



rameter t]s is unconstrained and we set it to - 1 for this fit only. 
Not surprisingly, most models are not able to match the lower- 
energy HEAO-3 data. For the lighter species, a + 6 remains the 
same, regardless of the model. The CREAM-II data for these 
species are > 100 GeV/n, in a regime where the asymptotic slope 
Tasympt IS reached. For the heavier species, where the data extend 
down to a few tens of GeV/n, a similar ordering (though less 
clear) of the models with yasympt (as for the HEAO-3 fit) is seen 
for the high-energy asymptotic behaviour. The fits to TRACER 
data shown in Fig. [TT] show an intermediate behaviour. Indeed, 
the energy range covers the same energy range as CREAM-II, 
but a few data point at low energy give a turn-over in the spec- 



trum. However, there is a gap between the two energy regimes, 
where the curvature of the HEAO-3 data is not reproducecfl. 

Goodness of fit Table|4]shows the best-fit spectral index a, and 
the associated Xmin/d o f ^alue for the various models, species, 
and data sets. Unsurprisingly, the best-fit (smaller A'^in/dof 
value) are for the CREAM-II data that only cover the high energy 



' Model I/O is not used on these data because it leads to unphysical 
values. This is likely to be related to the unphysical large negative value 
of r]s required, that becomes an issue for the numerical inversion of the 
energy losses (no convection and reacceleration to smooth out the steep 
upturn in the low-energy spectrum here, at variance with the three other 
models). 
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Table 4. Best- fit spectral index a and associated ^ ^ for the 
fit of the source spectrum parameters. Each column correspond 
to a different set of data on which the fit is performed. 



Model HEAO-3 TRACER CREAM 





d.o.f. 








Carbon 








TT 
il 


Z.4i 


7 1 fit 






9 zlS 


1 


in 


2.33 


6.14 


N/A 


N/A 


1.90 


2.15 


I/O 


2.28 


5.96 


N/A 


N/A 


2.11 


1.78 


111/11 


2.27 


6.54 


N/A 


N/A 


2.21 


1.73 








Oxygen 








TT 
1 1 


z..:> / 


7.34 


2.35 


1 1 7 
1 J. 1 / 


Z.OI 


^ 1 7 
J.LI 


TTT 


2.32 


6.11 


2.27 


13.76 


2.18 


4.59 


I/O 


2.26 


6.08 


N/A 


N/A 


2.28 


3.69 


111/11 


2.26 


6.61 


2.19 


2.54 


2.37 


3.54 








- Neon - 








TT 
1 1 


Z.J / 


3.94 


2.27 


J.J7 


9 S7 
Z.J / 


U.o J 


TTT 


2.30 


2.85 


2.19 


3.06 


2.01 


1.09 


I/O 


2.24 


3.54 


N/A 


N/A 


2.21 


0.93 


111/11 


2.23 


3.80 


2.09 


1.68 


2.31 


0.91 
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TT 
11 


Z.H-U 


7 1 T, 


Z,. J J 
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ZD. / 7 


9 SzL 
Z. J^- 


n AQ 


in 


2.35 


6.03 


2.23 


0.60 


2.10 


1.28 


I/O 


2.29 


6.37 


N/A 


N/A 


2.22 


0.90 


111/11 


2.29 


6.79 


2.17 


11.47 


2.31 


0.83 
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3.95 


2.49 


J J.OO 


9 An 
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9 ns 


m 


2.34 


3.15 


2.24 


3.89 


2.25 


2.93 


I/O 


2.29 


3.35 


N/A 


N/A 


2.33 


2.41 


111/11 


2.29 


3.63 


2.31 


35.14 


2.40 


2.30 
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9 99 
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VO 
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N/A 


N/A 


111/11 
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1.31 


2.19 


0.29 


N/A 
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2.56 


1.94 


2.49 


0.83 


N/A 


N/A 


m 


2.52 


1.76 


2.36 


1.08 


N/A 


N/A 


VO 


2.48 
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N/A 
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2.48 


1.90 


2.34 


0.73 


N/A 
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4.05 


2.48 


19.60 


2.67 


1.54 


ni 


2.39 


3.77 


2.29 


2.48 


2.42 


2.37 


I/O 


2.36 


3.88 


N/A 


N/A 


2.46 


1.94 


111/11 


2.35 


3.98 


2.34 


11.95 


2.52 


1.80 



range. It is indeed more difficult to reproduce the low energy 
part, where data have smaller error bars, but also where more 
effects (modulation, continuous and catastrophic losses) shape 
the spectrum. For the HEAO-3 case, the Xmm/d o f ^^1^^^ large 
for most of the species because of the difficulty to fit the highest 
energy point that has a very small error bar. For S, Ar and Ca, 
the fit is better. Data from the next CREAM flights, or from the 
AMS-02 instrument should help clarify the situation, and con- 
firm or otherwise these discrepancies amongst the various data. 
Nevertheless, some conclusions can still be drawn on the spec- 
tral indices (see below), although they are less constraining than 
those derived from the p and He data. 



Confidence intervals on a, q and rjs Figure [T2l shows, along 
with the 68% and 95% CIs, the best-fit values on the spectral 
indices of,- (top panel), the relative source abundances ^, (middle 
panel), and the source parameter 77^ , for all the primary species 
considered in this study. We first underline that the 95% CL rel- 
ative uncertainty for the parameters ranges from < 5% on a, and 
< 20% on <7,. CREAM data cover too naiTow an energy range to 
give stringent constraints, so we do not comment on them further 
below. The following trends are observed for the parameters: 

- a (top panel): as for the p and He data. Model II (reacceler- 
ation only, filled circle) always gives a larger value than the 
other models. Moreover, a similar range of slopes is found 
(2.2 - 2.5 for HEAO-3 data only, but more scatter when us- 
ing TRACER data). 

- qi (middle panel): the relative abundances from HEAO-3 
data are quite insensitive to the propagation model used. We 
recover the values of Engelmann et al. (199 0) (green boxes) , 
although with larger error bars. The lEngelmann et al.l(ll990h 
analysis is based on the leaky-box model, s o that their results 
are consistent with the iPutze et aP (l2009l) analysis (yellow 
boxes) performed in the same framework. Similar results are 
obtained in the diffusion model (this analysis), except for the 
discrepancy for S, Ar and Ca (our values are larger than those 
of the HEAO-3 analysis). The difference for Ar and Ca may 
be related to the fact these elements were assumed to be pure 
primary species in this analysis (in order to speed up the cal- 
culation). The lack of the secondary contribution translates 
in a higher primary flux required to match the data. This 
underlines the importance of taking properly into account 
all nuclei to derive the source abundances. Part of the dis- 
crepancy could also be related to the fact that, for the S, Ar, 
and Ca elements, the highest energy data point is better fitted 
than for the others (see above): this results in a larger value 
of rjs and a that may be responsible for the difference ob- 
served on the qj. The relative abundances obtained from the 
TRACER data are sensitive to the model chosen, presum- 
ably because of the lack of constraints in the intermediate 
energy range. The relative abundances for Model III (open 
squares) are consistent with those obtained from HEAO-3 
data, whereas the obtained values for Model II (filled cir- 
cles) and III/II (open diamonds) systematically undershoot 
those from HEAO-3 data. If we look into Tab. H) we remark 
that the fit to the Si data, on which all other abundances are 
normalised, is very poor for these models (large value of 
the , . t ). The consistently low values for all elements 
for these two models can be simply explained in terms of a 
too large Si abundance obtained from the TRACER data for 
Models II and III/II. The fact that Model III, which gives a 
good fit to Si data, gives abundances in agreement with those 
derived from HEAO-3 data, supports this explaination. 

- Tjs (bottom panel): as for the p and He data, a trend is ob- 
served showing a dependence on the models. The pattern is 
the same, but with the value of rjs larger by one than that for p 
and He. In terms of 775 -tjj, the C to Si primary species favour 
a value k 1, whereas p and He data favour k 0. There is more 
scatter in the S to Fe data, but we note that the HEAO-3 data 
are based on different use of sub-detectors for this heavier 
species. 

4.4. Summary for the source parameters 

An important result of the fixed transport parameter analysis is 
that, independent of the model and data considered, the source 
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Fig. 12. Best-fit value (symbols), 68% (dashed error bars) and 95% (solid error bars) CIs for C to Fe source parameters from the 
fit on CREAM-II, HEAO-3, and TRACER data . Top panel: source spectral index a,. Middle panel: relative abundances The 
results from the analys is of the HEAO-3 group dEngelmann e t al. 1990) are shown as green boxes, and those from the TRACER 
group (Ave et al. 2009) as orange boxes. The yellow boxes correspond to a leaky-box analysis of the abundances on HEAO-3 data 
performed in Putz e et al, (^00-9)- Bottom panel: low-energy source parameter rjs [see Eq.®]. Only the results from HEAO-3 data 
are plotted (for CREAM-II, rjs is set to -1, and for TRACER, the scatter is so large that the values are meaningless). 



slope of p and He nuclei is constrained to fall in the range 2.2 - 
2.5 (or 2.2-2.4 if we discard Model II). As the range of 6 covered 
by these models falls in the range 0.23 - 0.86 (see Table|2]i, this 
is a robust prediction. It also means that the asymptotic value for 



the propagated spectra (yasymp. = a + <5), that falls in the range 
2.7 - 3.0, is not reached in the GeV/n to TeV/n regime (as direct 
fits to the nuclei at Ek/n > 100 GeV/n propagated spectra lead 
to Tdata ~ 2.65. lAve et all I2OO8) : residual propagation effects 
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(reacceleration, convection, spallations) are still active. This is 
also supported by the results on the heavier species. 

Another important result is that regardless of the propagation 
model used, the quantity t/s - rjr is constrained to be ~ - 1 for 
all primary nuclei considered. Hence, to reproduce the data, if 
a pure power-law rigidity spectrum is assumed, a non-standard 
low-energy diffusion coefficient (upturn at low energy) is re- 
quired. Conversely, if a standard diffusion coefficient is assumed 
(i.e. K(E) = Kq/3'R.^), a flattening of the low-energy source spec- 
trum is required (i.e. dQ/dR oc with rjs > -1). The 
close-to-IS condition low-energy Voyager data is a further piece 
of information to break the degeneracy rjs - rfj, and would pos- 
sibly provide the shape of the low-energy source spectrum. 

Finally, most of the derived relative source abundances are 
in agreement with those derived by earlier groups. Still, a slight 
dependence on the propagation configuration is also observed. 
Due to the relevance of the value of the source abundances in 
the context of acceleration mechanisrns (.Dwver & Meveri,1987; 
IfVfever et alJ[l997t lEUison et alJll997t rOgliore et alJl2009l) . this 



point deserves further investigation. 

5. Ratio of primary species 

IWebber & LezniakI (11974 . more than 30 years ago, recognised 
the importance of looking at primary ratios. Such ratios may be, 
in principle, used to i) check the consistency of spectral indices 
of various species, ii) inspect whether source spectra are power- 
low in rigidity or power-law in kinetic energy, and also iii) in- 
spect whether solar modulation is a rigidity or total energy ef- 
fect. Below, we present several plots to illustrate some of these 
ideas, but also underline the complications that arise due to the 
many degeneracies between the source, transport, and modula- 
tion parameters (as underlined in the previous sections). 

5.1. p/He ratio 

Concerning the p/He raticQ on which IWebber & LezniakI (1 1974t) 
study mainly focused, the main conclusions were: i) proton and 
helium source spectra are rigidity rather than energy/nucleon 
spectra, ii) modulation effects dominate the shape of the p/He 
ratio for such rigidity spectra when shown as a function of ki- 
netic energy, and iii) modulation effects is not a pure rigidity 
effect since it flattens the spectrum at low energy relative to the 
interstellar flux. 



p/He from the toy-model calculation Caution is in order when 
calculating the ratio p/He, whether we start from the differen- 
tial flux in energy or in rigidity. In an analogous manner as for 
Eq. ©, 



dtj/ 
dEi 



k/n 



and 

_ # Z/3 
in order to define 



P_ 
He 



P_ 

He 



^ Note that because o f the misidentification of 'He and ''He, 
IWebber & LezniakI dl974) estimate a < 2.5% effect in the data plot- 
ting the same measurement as rigidity spectra or kinetic energy spectra, 
that is also not considered below. 



In the ID toy-model (energy gains and losses discarded), as- 
suming Eq. ^ for the source term — i.e. dQjdK - q/3''''^^'Rr" — , 
and Eq. (O for the diffusion coefficient — i.e. KCR) = 
P'^^ K{)R^ — , we have the analog of Eq. (|9l), but expressed in terms 
of the rigidity: 



P_ 
He 



2^p tii'^^m) 

/^l^jn) Ko'R'/{hL) + nc-/3l^J"{'R)-CTH. 
ySU' m Ko'R^/ihL) + nc ■ /sl'^K) ■ CTp 



(11) 



where we made explicit the rigidity dependence for all the terms 
(c is the speed of light). If the destruction rate is subdominant, 
we have 



He 



2qp /s';^'-'"('R) 



que 



(12) 



In all the above formulae, we have 



y8 = 



^jEklniEkjn + 2) 



For a proton jZ^ ^ 1, whereas ^ 4 for a helium nucleus. This 
is sufficient to distort the low energy p/He ratio whenever it is 
calculated from the differential fluxes in rigidity and rjs + 2 - 
TjT + 0. However, if the ratio is calculated from the differential 
fluxes in kinetic energy per nucleon, for a given Eki„, we have 
/^piEk/n) ~ P-aAEkin) = P, and Eq. (O reduces to 



P_ 

He 



^ 1p [Ek/„(Ek/„ + 2)] 



-(Q-p-aHe)/2 



X 



Ko[Ek/n(Ek/n + 2)]^/^2'V(/2L) + ncp'-''^cru, 
Ko[Ek/n(Ek/„ + iWl^lihL) + ncp'-nrcTp 

When the destruction rates are subdominant, we get 



He 



qp_ [Ek/„iEk/„+2)]-<-"^- "»'^/^ 



(13) 



(14) 



Comparison to data The p/He ratio is displayed as a function 
of the kinetic energy per nucleon in the left panel of Fig. [T3] 
for the BESS98 (red squares) and BESS-TeV (blue circles) data. 
The solid lines (no inelastic reaction terms) result from Eq. (fT4l i. 
with Q-p = aue- The shape of the ratio as well as the differences 
between the data taken at two different solar periods can be al- 
most completely ascribed to the modulation effect. The effect 
of the inelastic reaction term is contained in Eq. ( fTsT l. A closer 
look to this equation shows that the numerator and the denomi- 
nator do not differ by a factor of more ~ 3, confirming the sub- 
dominant (though important) role of this effect in determining 
the ratio (displayed versus kinetic energy per nucleon). The ef- 
fects (not shown) of having different spectral indices for p and 
He, having different values of a, 6, and rjs - rjj, when varied 
within reasonable limits, is of the same amount as the effect of 
the destruction rate. However, these effects are better seen when 
working with the rigidity. 

The right panel of Fig.[T3]shows a few experiments that have 
provided the p/He ratio as a function of the rigidity. Except for 
ATIC-1, the error bars were not provided by the experiments. 
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Fig. 13. Left panel: p/He ratio as a futiction of E^in, along with BESS98 (ISanuki et al.ll2000h and BESS-TeV (iHaino et alJl2004 
data. The lines show the toy-model calculation without the destruction term (thick solid lines) and with it (thick dotted lines). The 
black, red and blue lines are respectively modulated to O = MV (IS), O = 5 91 MV and O - 1 109. Left panel : same ratio, but 
as a function of the rigidity. The data are AMS-01 ( AMS Collaborati on et af] 12002). ATIC (Zatse pin etaP boOS). CAPRICE 94 
dBoezio et al.lll999i) . and some balloon data (Webber et al...l987) . The three toy-model calculations corresponds to the destruction 
rate set to zero (solid Unes) or to its required value (dotted Unes), plus a difference in the spectral index of p and He (dashed lines). 
In both plots, the normalisation is arbitrarily set to match the data. 



but are expected to be of the order of the size of the symbolfl 
As underlined in the previous sections, even though all experi- 
ments claim small error bars, not many of them are consistent 
with each other for the p and He fluxes. On the other hand, 
one would expect the ratio to have less systematics than fluxes. 
Nevertheless, a large discrepancy remains, which cannot be ex- 
plained by the different level of solar modulation associated to 
each experiment. The various curves show: i) the effect of solar 
modulation is sub-dominant for p/He vs H (thick vs thin lines); 
ii) the effect of inelastic interactions, which are switched off 
(solid lines) or included (dotted lines); iii) the distortion of the 
p/He ratio due to a possible difference in the spectral indices of p 
and He (dashed lines). The shape of the ratio depends mostly on 
the value of rjs - tjt- As found in Eq. (fTSI) . the ratio is constant 
if rjs - rjT = -2 (not shown on the figure). The best-fit to the 
data is obtained for rjs - tjt ~ I, in general agreement with the 
results of the more complete analysis of Sect. 14.21 A small dif- 
ference ftp - aue 5 0.05 between p and He source indices is not 
excluded. Note that the unperfect match to the data may result 
from the effect of energy losses which is not implemented in the 
toy formula. 

5.2. Ratio relative to tlie Oxygen fiux 

A similar analysis can be carried out for ratios of heavier species 
(Z > 2). In that case, for any element, we have A/Z ~ 2, so that 
for a given kinetic energy per nucleon, the rigidity or the /3 is the 
same for any element. The toy-model formulae is very similar to 
those for p/He Eq. ( fTST i: the parameters 775 and j]t, as well as the 

'The AMS-01 p and He flux given in l AMS Coflaboration et"al] 
(I2002h are not calculated for the same rigidity binning. We thus fit- 
ted both fluxes (expressed as a function of the nucleus rigidity) with 
a simple polynomial function, and calculated the ratio from these fits 
(shown as stars approximately along the same binning as the original 
He data^MS Collaboration et al. 2000). Calculating correctly the as- 
sociated error bars is not straightforward: as the data are just used for 
eye comparison, we do not go into more details. This difficulty in com- 
puting spectra binned in rigidity or kinetic energy suggests that in future 
experimenters should provide both. 



solar modulation effect are not expected to be important. This is 
shown in the top panel of Fig. [141 for a few elements: the main 
ingredient shaping the X/O ratios is the inelastic scattering on 
the ISM. 

There is a fair agreement with the data for C/0 (black), 
Si/O (orange) and Fe/O (magenta) as shown on the right 
panel of Fig. [141 especially at low-energy with the ACE data 
(iGeorge et all 12009) . The discrepancy with Ne/O and Mg/O is 
only at the level of ~ 20%. This could be due to systematics in 
the data, but this deserves further investigation, especially be- 
cause some isotopic anomalies in the Ne (and less likely for Mg) 
could be a signature for a contribution of the Wolf-Rayet stars 
to the standard cosmic-ray abundances ('Gupta & Webbeil l 19891 : 
Webber et al. 1997; Binns et al. 2005, 2008), or related to accel- 
eration in super-bubbles (e.g.. lHigdon & Lingenfelteij|2(j03h . 

6. Conclusions 

We have studied the spectra of proton, helium and other primary 
nuclei to derive the possible contraints and degeneracies on the 
source spectrum and the transport parameters. We have checked 
the compatibility of the primary fluxes with the transport pa- 
rameters derived from the B/C analysis, and inspected whether 
they add further constraints. We have then derived the source pa- 
rameters: slope of the power-law spectrum, abundance and low- 
energy shape. 

In Sect. 3, we have analysed the fluxes of primary cosmic 
rays in diffusion models with particular attention to p and He. 
The most recent data on p and He are well reproduced by a 
purely diffusive model, described by power-law source spec- 
trum, isotropic diffusion coefficient, spallative destructions and 
electromagnetic energy losses. This conclusion holds for single 
data sets but it is not reproduced in a combined data analysis 
(except for AMSOl and BESS98 proton data), due to the mean 
level of consistency among the different data collections. The in- 
spection of low energy (< 100 GeV/n) p and He data indicates 
that the purely diffusive regime is likely not reached due to the 
role of spallations and, to less extent, of energy losses. The data 
are shown to be compatible with a wide class of purely diffusive 
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Fig. 14. Ratio of element to O as a function of the kinetic en- 
ergy per nucleon. Top panel: C/O, Si/O and Fe/O without (IS, 
solid lines) or with (TOA, dashed lines) solar modulation. The 
destruction cross-section is indicated for each element. Bottom 
panel: comparison of the simple toy-model for mula (thick solid 
lines) with the data (symbols): CREAM I and II ("Ahn et al."2008', 
[2009, 2010a b), CRNjSwordv et al. 1990; Mueller et al. 1991) 
HEAO -3 dEngehnann et all Il990t). TRACER ([Ave et all l2008i 
120091) . and low-energy ACE data dGeorge et al.ll2009h . 



models, but can also be accommodated by models with convec- 
tion and/or acceleration. In all scenarios, they do not put signif- 
icant constraints on the transport parameters and tend to favour 
values for the source and transport parameters outside reason- 
able physical limits. 

In Sect. 4 we have applied best-fit models previously selected 
from B/C data, to the propagation of light and heavier primary 
nuclei, where we consider the possibility of low-energy devia- 
tions both in the diffusion coefficient K{K) - KqP^'^'R? and in 
the acceleration spectrum dQ/dR = qfS'^^^^'R'" . The main re- 
sults on p and He are: i) it is possible to accommodate these pri- 
mary fluxes in diffusive models along with B/C data (very good 
fit on BESS data, less satisfactory on AMS ones); ii) a ranges 
between 2.2 and 2.5; iii) for any given data set and species the 
spread in the source power index is Up - q-hc S 0.2, regardless 
of the model; iv) p and He point to very similar rjs , whose val- 
ues depend on the model: close to 1 for the reacceleration and 
convection/reacceleration models, whereas rjs is close to -2 for 
Model I/O and -1.5 for Model III/II, which contain a low energy 
upturn in the diffusion coefficient by means of 7/7. Indeed, the 
constraint which seems to emerge from any propagation model 



is on their difference: 77s - //j- w - 1. We have demonstrated 
that a possible way to break the degeneracy rjs - t]t inherent the 
low-energy tail is by means of the Voyager data, taken at a few 
hundreds of MeV/n, and in a quasi-IS regime. 

We have also studied heavier primary nuclei spectra, whose 
relevant destruction rate on the ISM increases roughly with 
atomic number. Therefore, they can be tested against a univer- 
sal source spectral function. As for p and He, we have fitted data 
for C, O, Ne, Mg, Si, S, Ar, Ca and Fe using the models selected 
based on B/C. As for light primaries, these nuclei point to a 
around 2.2-2.5 regardless of the propagation configuration used. 
The quantity rjs - rfj is also constrained from C to Fe primary 
data to be ?7s - 777- ^ 1 . Moreover, most of the derived relative 
source abundances are in agreement with those derived by ear- 
lier groups. We have shown that it is possible to reproduce the 
primary cosmic -ray spectra from protons to iron for all consid- 
ered models (with or without convection/reacceleration) without 
a need for an artifici al break in the injec tion slope at 10 GV, as 
it is assumed, e.g., in lTrotta et al. ( 2OI0h. This is also consistent 
with the leaky-box analysis of Putze et aP (120091 

In Sect. 5, we have studied the ratio of two primary species. 
The shape of the p/He ratios, when plotted as a function of the 
kinetic energy per nucleon, is driven by the modulation effect. 
This effect could be used to monitor the modulation level at dif- 
ferent periods. On the other hand, the same ratio plotted as a 
function of the rigidity minimises the effect of the modulation, 
and is well adapted to probe the values of rjs - tjt, and ap - aue- 
A full analysis using energy gains and losses is required to fur- 
ther develop such an approach, which is complementary to the 
direct fit of the p and He fluxes, and which may suffer less from 
systematics in the data. For X/O ratios, where X is an primary el- 
ement with Z > 2, the behaviour as a function of Ek/n is driven 
by the destruction rate on the ISM. Isotopic anomalies and/or 
non-universality of the source slopes can be inspected by means 
of these ratios. 

For the primary spectra and ratio of primary spectra, the 
importance of the destructive terms on the measured spectra 
has not previously been much considered. This is an innova- 
tion of the paper. We have also detailed the degeneracy between 
the source and propagation parameters, and provided probabil- 
ity density functions for the spectral index, abundance and low- 
energy shape of the source spectra for all species, as a means to 
address the question of the universality of the source spectra and 
their low-energy shape. Our analysis reinforces the need of more 
accurate data on light primary nuclei not only in the low-energy 
regime but also in the TeV/n-PeV/n range, as well as accu- 
rate measurement of primary-to-primary ratio. With such data 
it would be possible to significantly constrain the low-energy 
shape of the diffusion coefficient and the source spectrum, and 
fix the asymptotic behaviour of propagated nuclei. 
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Appendix A: Source slopes and relative elemental abundances 



Table A.l. Most probable values and 68% CIs on the spectral index a and the relative source abundance qzl<lsi for the fit of the 
source spectrum parameters. Each column corresponds to a different set of data on which the fit is performed. 
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